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We report the first observation of the spin-singlet bottomonium states h b (lP) and h b (2P) produced 
in the reaction e + e~ — > h b {nP)ix + -K~ using a 121.4 fb _1 data sample collected at energies near the 
T(5S) resonance with the Belle detector at the KEKB asymmetric-energy e + e~ collider. We deter- 
mine M\h b (lP)] = (9898.3 ± 1.1^;?) MeV/c 2 and M[h b (2P)] = (10259.8 ± O.Gtl f,) MeV/c 2 , which 
correspond to P-wave hyperfine splittings AM HF = (+1.6 ± 1.5) MeV/c 2 and (+0.5±i;|) MeV/c 2 , 
respectively. The h b (lP) and h b (2P) are observed with significances of 5.5 a and 11.2 a, respec- 
tively. We also report measurements of the cross sections for e + e~ — > h b (nP)-K + -K~ relative to that 
for e+e- -> T(2S)tt + tv'. 

PACS numbers: 14.40. Pq, 13.25. Gv, 12.39.Pn 



Bottomonium is the bound system of bb quarks and 
is considered an excellent laboratory to study Quantum 
Chromodynamics (QCD) at low energies. The system is 
approximately non-relativistic due to the large b quark 
mass, and therefore the quark-antiquark QCD potential 
can be investigated via bb spectroscopy 

The spin-singlet states h b (nP) and rj b (nS) alone pro- 
vide information concerning the spin-spin (or hyper- 
fine) interaction in bottomonium. Measurements of the 
h b (nP) masses provide unique access to the P-wave hy- 
perfine splitting, AM H f = (M(n 3 Pj)) - Af(n 1 Pi), the 
difference between the spin-weighted average mass of the 
P-wave triplet states (xbj(nP) or n 3 Pj) and that of the 
corresponding h b (nP), or n 1 Pi. These splittings are pre- 
dicted to be close to zero 0| , and recent measurements of 
the h c (lP) mass correspond to a P-wave hyperfine split- 
ting that validates this expectation for the IP level in 
charmonium: AM HF = (0.00 ± 0.15) McV/c 2 Q. 

Recently, the CLEO Collaboration observed the pro- 
cess e + e~ — > h c (lP)Tr + ir~ at a rate comparable to that 
for e + e~ — > J/ipTr + ir~ in data taken above open charm 
threshold Such a large rate was unexpected be- 

cause the production of h c {lP) requires a c-quark spin- 
flip, while production of J ftp does not. Similarly, the 
Belle Collaboration observed anomalously high rates for 
e + e~ — > T(nS)Tr + iv~ (n = 1,2,3) at energies near the 
T(55) mass Q. Together, these observations motivate 



a search for e + e — > n + ir h b (nP) above open-bottom 
threshold at the T(55) resonance. 

In this Letter, we report the first observation of the 
h b (lP) and h b (2P) produced via e + e~ — > h b (nP)ir + ir~ 
in the T(55*) region. We use a 121.4 fb _1 data sample 
collected near the peak of the T(55) resonance (y/s ~ 
10.865 GeV) with the Belle detector @ at the KEKB 
asymmetric-energy e + e _ collider Q. 

We observe the h b (nP) states in the 7r + 7r~ missing 
mass spectrum of hadronic events. The 7r + 7r _ missing 
mass is defined as M^ liss = (Pr(5S) — Pn+K-) 2 , where 
Pt(5S) is the 4-momentum of the T(5S) determined from 
the beam momenta and P^^- is the 4-momentum of 
the 7r + 7r _ system. The 7r + 7r~ transitions between T(nS) 
states provide high-statistics reference signals. 

Our hadronic event selection requires a reconstructed 
primary vertex consistent with the run-averaged inter- 
action point (IP), at least three high-quality charged 
tracks, a total visible energy greater than 0.2 y/s, a total 
neutral energy of (0.1 — 0.8) y/s, more than one large- 
angle cluster in the electromagnetic calorimeter and that 
the total center-of-mass momentum have longitudinal 
component smaller than 0.5 yfs [8j. The 7r + 7r~ candi- 
dates are pairs of well reconstructed, oppositely charged 
tracks that are identified as pions and do not satisfy 
electron-identification criteria. Continuum e + e~ — > qq 
(q = u, d, s, c) background is suppressed by requir- 
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FIG. 1. The Af m i ss distribution for the selected 7r + 7r~ pairs. 
Vertical lines indicate the locations of the T(IS'), hb(lP), 
T(2S), h b (2P) and T(3S) signals. 

ing the ratio of the second to zeroth Fox- Wolfram mo- 
ments to satisfy Ri < 0.3 The resulting M m i ss spec- 
trum, which is dominated by combinatoric 7r + 7r~ pairs, 
is shown in Fig. Q] 

Prior to fitting the inclusive M m i SS spectrum we study 
reference channels and peaking backgrounds arising from 
7r + 7r~ transitions between T(nS) states. A high purity 
sample of such transitions is obtained by reconstructing 
pairs in the event in addition to the n + ir~ pair. 
For these studies the hadronic event selection criteria are 
not applied, while for the /x + fi~ pair we use the same 
selection as was employed in Ref. @. MC studies indi- 
cate that the shape of the peaks in M m i ss is independent 
of whether the tt + tt~ arc reconstructed in the hadronic 
environment or in this much cleaner environment. In ad- 
dition, to suppress radiative Bhabha events in which the 
photon converts, producing a fake 7r + 7r~, we require that 
the opening angle between the candidate pions in the lab- 
oratory frame satisfies cos0 ff + ff - < 0.95. In Fig. [5] (a) we 
present the two-dimensional distribution of mass 
A^u+u- vs - -^miss for events satisfying these criteria. 

Clear peaks are visible along a diagonal band, where 
M^+p- is roughly equal to -M m i ss , and correspond to 
fully reconstructed T(55) — !• T(nS)n + Tr~ — > /i + ^~7r + 7r~ 
events. Also along the diagonal is a diffuse background of 
events that arise due to the process e + e" — > 
e + e _ ), where the conversion pair is reconstructed as 
7r + 7r~, or from non-resonant e + e~ — > fi + fi~it + n~ events. 
Events from the band satisfying |M m ; ss — M (Jt + (Jt -\ < 
150MeV/c 2 are projected onto the M m ; ss axis and fit- 
ted to the sum of a linear background and a Gaussian 
joined to a power-law tail on the high mass side. The 
high-side tail is due to Initial State Radiation (ISR) pho- 
tons. This latter function is analogous to the well-known 
Crystal Ball function but has the tail on the higher 
rather than lower side. We thus refer to it as a 'reversed 
Crystal Ball' (rCB) function. The fitted M m i ss spectra 



TABLE I. The yield, mass and width for signals recon- 
structed using M m i S s from the exclusive p + fi~ tt + tt~ selection. 
Each mass is consistent with the world average [ill ]. 





Yield 


Mass, MeV/c 2 


a, MeV/c 2 


T(15) 


1894 ± 61 


9459.96 ± 0.23 


7.68 ± 0.21 


T(25) 


2322 ± 60 


10023.34 ± 0.22 


6.60 ± 0.20 


T(3S) 


661± 3 3 9 


10355.66±g;|| 


5.98±g;«? 



from this band are shown in Figs. [H (b)-(d), and the re- 
sulting yields, masses and width of the rCB function for 
the T(nS) states are displayed in Tabled] The masses ob- 
tained are consistent with the world average values [ll| . 

The structures in the horizontal band in Fig. [5] (a), 
where M M + M - is roughly equal to M[T(15)], arise from 
events in which a daughter T(IS') in the event de- 
cays to In Figs. [2] (e)-(f) we present M m j ss 
projections from this band, subject to the requirement 
- M[T(IS)]\ < 150MeV/c 2 . The peaks at the 
T(35) and T(25) masses arise from events having 7r + 7r~ 
transitions to T(3S) or Y(25), followed by inclusive pro- 
duction of T(15), and are fitted to rCB functions. Peaks 
at 9.97GeV/c 2 and 10.30 GeV/c 2 arise from events in 
which a T(35) or T(25) is produced inclusively in T(55) 
decays or via ISR, and then decays to T(lS)ir + n~ , and 
are fitted to single and double Gaussians, respectively. 

The threshold for inclusive Kg production results in a 
sharp rise in the M m i ss spectrum, due to Kg — > tt + tt~, 
very close to the mass of T(35). Rather than veto ir + Tr~ 
combinations with invariant masses near M(Kg), which 
significantly distorts the M m i SS spectrum in the vicinity, 
we obtain the Kg contamination by fitting the ir + 7r~ 
invariant mass corresponding to bins of Af m ; ss . 

The M m ; ss spectrum is divided into three adjacent re- 
gions with boundaries at Af m i SS = 9.3, 9.8, 10.1 and 
10.45 GeV/c 2 and fitted separately in each region. In 
the first two regions, we use a 6th-order Chebyshev poly- 
nomial, while in the third we use a 7th-order one. In the 
third region, prior to fitting, we subtract the contribution 
due to Kg — > tt + tt~ bin-by-bin. The signal component 
of the fit includes all signals seen in the (JT tt + tt~ data 
as well as those arising from tt + it~ transitions to hb(nP) 
and T(1D). We fit these additional signals using the tail 
parameters of the T(2S) and fixed widths found by linear 
interpolation in mass from the widths of the exclusively- 
reconstructed T(nS) peaks. The peak positions of all sig- 
nals are floated, except that for T(3S I ) — > Y(1S)tt + it~ , 
which is poorly constrained by the fit. The confidence 
levels of the fits in the three regions are 3.0%, 0.5% and 
0.4%, respectively. The M m ; ss spectrum, after subtrac- 
tion of both the combinatoric and Kg — > tt + it~ contri- 
butions is shown with the fitted signal functions overlaid 
in Fig. [31 The signal parameters are listed in Tabic HT1 

We studied several sources of systematic uncertainty. 
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FIG. 2. (a) Distribution of Af +„- vs. M m i ss , and the projection on M m i aa from (b)-(d), the diagonal band |Af m i ss — M„+„- \ < 



150MeV/c 2 near the T(1S), T(2S) and T(3S); and (e),(f), the horizontal band \M^ 
T(2S) and T(35). 



M[T(15)]| < 150MeV/c 2 near the 



TABLE II. The yield, mass and statistical significance from 
the fits to the M m i ss distributions. The statistical significance 
is calculated from the difference in x 2 between the best fit and 
the fit with the signal yield fixed to zero. 





Yield, 10 3 


Mass, MeV/c 2 


Significance 


T(1S) 


105.2 ±5.8 ±3.0 


9459.4 ± 0.5 ± 1.0 


18.2 a 


h b {lP) 


50.4 ± 7.8±_ti 


9898.3 ±1.111;? 


6.2 a 


3S^ IS 


56 ±19 


9973.01 


2.9 a 


T(2S) 


143.5 ±8.7 ±6.8 


10022.3 ±0.4 ±1.0 


16.6(7 


T(1L>) 


22.0 ± 7.8 


10166.2 ± 2.6 


2.4 a 


h b (2P) 


84.4 ± 6.8tlo. 


10259.8 ± Q.6±{-% 


12.4 a 


2S -> IS 


151.7 ± 9.7^o 


10304.6 ± 0.6 ± 1.0 


15.7a 


T(3S) 


45.6 ±5.2 ±5.1 


10356.7 ± 0.9 ± 1.1 


8.5 a 



The background polynomial order was increased by 
three, and the range of the fits performed were altered by 
up to 100 MeV/c 2 . Different signal functions were used, 
including symmetric Gaussians and rCB functions with 
the width parameters left free. We altered our selection 
criteria: tightening the requirements on the proximity of 
track origin to the IP, increasing the minimum number 
of tracks to four, and imposing the cos^+^r- < 0.95 re- 
quirement used in the /i + /i~7r + 7r~ study. In Table ITTT1 a 
summary of our systematic studies is presented. 

The values in the table represent the maximal change 
of parameters under the variations explored. We esti- 
mate an additional 1 MeV/c 2 uncertainty in mass mea- 



TABLE III. Absolute systematic uncertainties in the yields 
and masses from various sources. 





Polynomial 
order 


Fit 
range 


Signal Selection 
shape requirements 


A[T(15')], 10 3 


±1.4 


±1.7 


±2.0 




N[h b (lP)], 10 3 


±2.4 


±3.6 


+ 1.2 
-8.0 




A[T(2S')], 10 3 


±3.4 


±3.2 


±5.0 




N[h b (2P)], 10 3 


±2.2 


±2.6 


+23. 
-9.0 




N[2 -> 1], 10 3 


±3.0 


±8.0 


+0 
-18 




AT[T(35)], 10 3 


±1.0 


±3.0 


±4.0 




M[T(15")], MeV/c 2 


±0.04 


±0.06 


±0.03 


±0.18 


M{h b (lP)], MeV/c 2 


±0.04 


±0.10 


+0.04 
-0.20 


+0.20 
-0.30 


M[T(2S)], MeV/c 2 


±0.02 


±0.08 


±0.06 


±0.03 


M[h b (2P)}, MeV/c 2 


±0.10 


±0.20 


+ 1.0 
-0.0 


±0.08 


M[2 ->■ 1], MeV/c 2 


±0.20 


±0.10 


±0.06 


±0.10 


M[T(3S)], MeV/c 2 


±0.15 


±0.24 


±0.10 


±0.20 



surements based on the differences between the observed 
values of the fitted T(nS) peak positions and their world 
averages. The total systematic uncertainties presented in 
Table [TT] represent the sum in quadrature of all the con- 
tributions listed in Table UIT1 The signal for the T(1Z?) 
is marginal and therefore systematic uncertainties on its 
related measurements are not listed in the table. The 
significances of the /ib(lP) and h\,{2P) signals, with sys- 
tematic uncertainties accounted for, are 5.5(7 and 11.2cr, 
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FIG. 3. The inclusive A/ m i ss spectrum with the combinatoric background and K% contribution subtracted (points with errors) 
and signal component of the fit function overlaid (smooth curve). The vertical lines indicate boundaries of the fit regions. 



respectively. 

The measured 
M = (9898.3 ± 1.1 



masses 

1.0 

i.i 



of h b (IP) and 
MeV/c 2 and M = 



h b (2P) are 
(10259.8 ± 



0.6_ 1 ' ) MeV/c , respectively. Using the world average 
masses of the Xbj(nP) states, we determine the hyper- 
fine splittings to be AAf H F = (+1-6 ± 1.5)MeV/c 2 and 
(+0.51^2) MeV/c 2 , respectively, where statistical and 
systematic uncertainties are combined in quadrature. 



We also measure the ratio of cross sections for 



T(5S) -> h b (nP)TT+TT- to that for e+er -> T(5S) -> 
T(2S)tt + it~ . To determine the reconstruction effi- 
ciency we use the results of resonant structure stud- 
ies reported in Ref. 



121 that revealed the existence 



of two charged bottmonium-like states, Zb(10610) and 
•Z&(10650), through which the 7r + 7r~ transitions we are 
studying primarily proceed. These studies indicate that 
the Z b most likely have J p = 1 + , and therefore in our 
simulations the ir + ir~ transitions are generated accord- 
ingly. To estimate the systematic uncertainty in our re- 
construction efficiencies, we use MC samples generated 
with all allowed quantum numbers with J < 2. 

We find that the reconstruction efficiency for the 
T(2S) is about 57%, and that those for the h b (lP) and 
h b (2P) relative to that for the T(25) are 0.913^ n?n and 



0.824 



+0.130 
-0.013' 



0.010 

respectively. The efficiency of the R2 < 0.3 



requirement is estimated from data by measuring signal 
yields with R 2 > 0.3. For T(2S), h b (lP) and h b (2P) 
we find 0.863 ± 0.032, 0.723 ± 0.068 and 0.796 ± 0.043, 
respectively. From the yields and efficiencies described 
above, we determine the ratio of cross sections R = 
$7Wj£fj to be i? = 0.46 ± 0.08l°;?2 6» the h b (lP) 
and R = 0.77±0.08^; 22 for the h b (2P). Hence T(55) 
h b (nP)ir + ir~ and T(55) — > T(2S')7r + 7r _ proceed at sim- 



ilar rates, despite the fact that the production of h b (nP) 
requires a spin-flip of a &-quark. 

The rate of T(55) — > h b (nP)ir + ir~ is much larger than 
the upper limit for that of T(3S I ) — > h b (nP)ir + ir~ ob- 
tained by the BaBarCollaboration [ljj . This is consis- 
tent with the observation that the rates for T(55) — > 
Y (mS)Tt + n~ with m = 1,2,3 are much larger than those 
for T(nS) -> r{mS)ir+ir- for n = 2,3,4 @. The 
only previous evidence for the h b (lP) is a 3. Oct excess 
in T(35) -> n°h b (lP) at (9902 ± 4) MeV/c 2 presented 
by BaBar [3]. 

We have also used 711 fb — 1 of e + e~ collisions at the 
T(4.S) resonance to search for T(45) — > h b (lP)ir + ir~ 
(h b (2P) is kinematically forbidden). The overall effi- 
ciency, assuming the R2 efficiency at T(45) to be the 
same as that at T(55), is 0.94jlQ Q3 relative to that for 
T(55) — > h b (lP)-K + iT~ . From our observed yield of 
(35 ± 2lt\t) x 10 3 , we therefore set an upper limit on 
the ratio of CT(e+e~ -> h b (lP)n+TT-) at the T(45) to 
that at the T(5S') of 0.27 at 90% C.L. 

In summary, we have observed the P-wave spin-singlet 
bottomonium states h b (lP) and h b (2P) in the reaction 
e+er -> T(55) -> h b (nP)ir + ir- . The h b {nP) masses 
correspond to hyperfine splittings that are consistent 
with zero. We also have observed that the cross sec- 
tions for these processes and that for e + e~ — > T(5S) — > 
T(2S f )7r + 7r~ are of comparable magnitude, indicating the 
production of h b (nP) at the T(5S) resonance must oc- 
cur via a process that avoids the expected suppression 
related to heavy quark spin-flip. 
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